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Abstract 

We consider a three-dimensional chaotic system consisting of the suspension of Arnold’s 
cat map coupled with a clock via a weak dissipative interaction. We show that the coupled 
system displays a synchronization phenomenon, in the sense that the relative phase between 
the suspension flow and the clock locks to a special value, thus making the motion fall onto a 
lower dimensional attractor. More specifically, we construct the attractive invariant manifold, 
of dimension smaller than three, using a convergent perturbative expansion. Moreover, we 
compute via convergent series the Lyapunov exponents, including notably the central one. The 
result generalizes a previous construction of the attractive invariant manifold in a similar but 
simpler model. The main novelty of the current construction relies in the computation of the 
Lyapunov spectrum, which consists of non-trivial analytic exponents. Some conjectures about a 
possible smoothening transition of the attractor as the coupling is increased are also discussed. 

Keywords: partially hyperbolic systems; Anosov systems; synchronization; phase-locking; Lya¬ 
punov exponents; fractal attractor; SRB measure; tree expansion; perturbation theory. 


1 Introduction 

Synchronization in chaotic systems is a surprising phenomenon, which recently received a lot of 
attention, see e.g. Even though the heuristic theory and the classification of 

the synchronization phenomena are well studied and reasonably well understood, a mathematically 
rigorous theory is still lacking. Generally speaking, a standard difficulty lies in the fact that the 
phenomenon involves the dynamics of non-uniformly chaotic systems, typically consisting of different 
sub-systems, whose long-time behavior depends crucially on the sign of the “central” Lyapunov 
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exponents, i.e. of those exponents that are zero in the case of zero coupling, and become possibly 
non-trivial in the presence of interactions among the sub-systems. The mathematical control of such 
exponents is typically very hard. Progress in their computation is a fundamental preliminary step 
for the construction of the SRB measure of chains or lattices of chaotic flows, which may serve as 
toy models for extensive chaotic systems out-of-equilibrium (i.e. they may serve as standard models 
for non-equilibrium steady states in non-equilibrium statistical mechanics). 

In a previous paper [9], we introduced a simple model for phase synchronization in a three- 
dimensional system consisting of the suspension flow of Arnold’s cat map coupled with a clock. 
The coupling in [9] was unidirectional, in the sense that it did not modify the suspension flow, but 
only the clock motion. Notwithstanding its simplicity, the model has a non-trivial behavior: in 
particular, it exhibits phase locking and in [9j we constructed the corresponding attractive invariant 
manifold via a convergent expansion. However, because of unidirectionality, the Lyapunov spectrum 
in [9] was very simple: the “longitudinal” exponents (i.e., those corresponding to the motion on the 
invariant manifold) coincided with the unperturbed ones, and the central exponent was expressed 
in the form of a simple integral of the perturbation over the manifold. In this paper, we extend the 
analysis of |9j to a simple bidirectional model, for which the Lyapunov spectrum is non-trivial, and 
we show how to compute it in terms of a modified expansion, which takes the form of a decorated 
tree expansion discussed in detail in the following. 

The model is defined as follows. Take Arnold’s cat map x i —> Sx and denote by A-t and x± the 
eigenvalues and eigenvectors, respectively, of S : 


S = 




A± = 


3 ± \/5 



( 1 . 1 ) 


with A f± := \Jl + (A± — 2) 1 2 , so that x± are normalized. We let the suspension flow of Arnold’s cat 
be defined as x —> x{t) := d> f (x), with x{t) = S k x, if 2irk <t< 2n(k + 1) V7c € Z. Formally, x(t) is 
the solution to the following differential equatioi^j] on T 2 = M 2 /(27rZ): 


x = 5(t) (log S)x , (1.2) 

where 5 is the 2-7r-periodic delta function such that J 0 * dt'5(t') = 1 for all 0 < t < 27r. The model of 
interest is obtained by coupling the suspension flow of Arnold’s cat map with a clock by a regular 

1 In 0 we erroneously wrote x = 8(t){S — l)x instead of B but throughout the paper we only used the fact 
that at all times 2nk + 0 + the variable x jumped abruptly from x(2nk) to Sx(2nk), and besides these discontinuities 
the flow was smooth. Therefore, all the results and statements of [5] are correct, modulo this re-interpretation of the 

flow equation [9] (2.1)], where S(t)(Sx — x) should be replaced by 8(t)(log S)x. 
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perturbation, so that on T 2 X T the evolution equation is 


x = 5(t) [ log S~\x + sf(x, w, t), 
w = 1 + eg(x, w, t ), 


(1.3) 


where e > 0 and /, g are 27r-periodic in their arguments. For e = 0 the motions of x and w are 
independent. Therefore, the relative phase w mod 2ir among the two flows is arbitrary. If e 0 and 
if the interaction is dissipative (in a suitable sense, to be clarified in a moment), then the phases 
of the two sub-systems can lock, so that the limiting motion in the far future takes place on an 
attractor of dimension smaller than 3, for all initial data in an open neighborood of the attractor. 
In [9], we explicitly constructed such an attractor in terms of a convergent power series expansion 
in e, for / = 0 and a special class of dissipative functions g. In this paper, we generalize the analysis 
of [9] to / 0. Our first result concerns the construction of the attractive invariant manifold for 

ffl o. 

Theorem 1.1. Let S * be the flow on T 3 associated with the dynamics 
in their arguments. Set 

r2n r2n 

1 q(<P,w q ):= dtg(Sip,w 0 + t,t), 71 (<p, w 0 ) := d td w g(S<p,w Q + t,t), (1.4) 

Jo Jo 

and assume there exists wo G T such that 70 ( 92 , wq) = 0 and T := 71(97 tco) < 0, independently 
of ip. Then there are constants £q, c > 0 such that for 0 < e < Eq there exist a homemorphism 
H : T 2 —> T 2 and a continuous function W : T 2 — > T, both Holder-continuous of exponent fl > ce, 
such that the surface II = {(H(ip),W(tp)) : ip G T 2 } is invariant under the Poincare map and 
the dynamics of on f 1 is conjugated to that of S on T 2 , i.e. 


(1.3), with f and g analytic 


S 2 /(H(<p),W(<p)) = (H(Sip),W(Sip)). 


(1.5) 


The proof of this theorem is constructive: it provides an explicit algorithm for computing the 
generic term of the perturbation series of (H, W ) with respect to e, it shows how to estimate it and 
how to prove convergence of the series. As a by-product, we show that the invariant manifold is 
holomorphic in £ in a suitable domain of the complex plane, whose boundary contains the origin. 
The construction also implies that Ll is an attractor. We denote by B its basin of attraction and by 
Ho an arbitrary open neighborood of Ll contained in B such that plq{Bq) > 0, with (j,q the Lesbegue 
measure on T 3 . 

In addition to the construction of the invariant surface, in this paper we show how to com- 
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pute the invariant measure on the attractor and the Lyapunov spectrum, in terms of convergent 
expansions. More precisely, let p be the Lesbegue measure restricted to Bq, i.e., denoting by xb 0 
the characteristic function of Bq, fi(A) := po(xb 0 A), for all measurable A. The “natural” invariant 
measure on the attractor, /JSRB■ is defined by 

lim j f F(Sl{x))dT = iiskb(F) 

t^oo t J o 

for all continuous functions F and p-a.e. x, where x = (x,w) 6 T 3 . The limiting measure //srb is 
supported on and such that //srb(^) := /^SRB(xn) = 1- On the attractor, /USRB-a.e. point defines 
a dynamical base, i.e. a decomposition of the tangent plane as T T 3 = W+(x) © Wo(x) © W_(x), 
such that 

1 |<9S 27 rfe (xW£| 

A lirn ^ log ---- = A© if 0 / d£ <E Wj(x) , i G {+, 0, -} . 

The constants of motion Aj are the Lyapunov exponents, and we suppose them ordered as A + > 
Ao > A_; in the following we shall call Ao the central Lyapunov exponent. Our second main result 
is the following. 

Theorem 1.2. There exists eq such that the following is true. Let F be an Holder continuous 
function on T 3 . Then pskb (F) is Holder continuous in e, for 0 < e < £q. If F is analytic, then 
hSRB(F) is analytic in e, for 0 <£<£f<£q and a suitable F-dependent constant £p. Moreover, 
the Lyapunov exponents A,, i e {+, 0, —} ; are analytic in £ for 0 < £ < £q. In particular, the central 
Lyapunov exponent is negative: Ao = eT + 0(e 2 ), while A± = log X± + 0(e). 

The paper is organized as follows. Theorem 0 is proved in Section [2] below. The proof follows 
the same strategy of [9]: (1) we first write the equations for the invariant surface and solve them 
recursively at all orders in e; ( 2 ) then we express the result of the recursion (which is not simply 
a power series in e) in terms of tree diagrams (planar graphs without loops); trees with n nodes 
are proportional to £ n times a tree value, which is also a function of e; (3) finally, using the tree 
representation, we derive an upper bound on the tree values. The fact that the dissipation is small, 
of order 0 (e), produces bad factors E~ k in the bounds of the tree values, for some k depending on 
the tree. Therefore, we need to show that for any tree k is smaller than a fraction of n, if n is the 
number of nodes in the tree. This is proved by exhibiting suitable cancellations, arising from the 
condition 70 (</?, wo) = 0 . 

Theorem |1.2| is proved in Section [3] The proof adapts the tree expansion to the computation of 
the local Lyapunov exponents Aj(x) on the invariant surface, in the spirit of [5! Chapter 10] and 0. 
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The positive local Lyapunov exponent A+(x) plays the role of the Gibbs potential for the invariant 
measure g>sRB- Therefore, given a convergent expansion for A+(x), p srb can be constructed by 
standard cluster expansion methods, as in Chapter 10]. Finally, A, can be expressed as averages 
of the local exponents over the stationary distribution. 

In Section [4j we present some numerical evidences for a fractal to non-fractal transition of the 
invariant manifold, and formulate some conjectures. 


2 Formation of an invariant surface 

2.1 Conjugation 

In this section, we define the equations for the invariant manifold, by introducing a conjugation 
that maps the dynamics restricted to the attractor onto the unperturbed one. The conjugation is 
denoted by (H, W) : T 2 —» (T 2 ,T), with 


H(ip) :=ip + h((p) = (l + /i)(^), W(tp) :=w 0 + U(ip), 


where 1 is the identity in M 2 . 

Let x(0) = x = H(ip) and w( 0) = wq + U(<p) be the initial conditions at time t = 0. We will 


look for a solution to (1.3) of the form 


x(t) = Sx + a(ip, t ) = SH(ip) + a(ip, t ) = Sip + Sh(ip) + a(ip, t), 
w(t) = w 0 + t + 


( 2 . 1 ) 


for 0 < t < 2n, with boundary conditions 


a((/?,0 + ) = 0, a(ip, 2tt) = H(Sip) — SH(tp), u(ip, 0 + ) = U(ip), u(ip, 2n) = U(S<p). (2.2) 


The evolution equation for w will be written by “expanding the vector field at first order in u and 
at zeroth order in Sh + a”, i.e. as 

u(ip,t) = eg(x(t),w{t),t) = £7o (<P,t) +£71 (ip,t) u{<p,t) +eG(<p,t), (2.3) 

where 

:= g{S(p,w 0 + t,t), 7i (<p,t) := d w g(Stp,w 0 + t,t) (2.4) 
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and 


G(<p,t) := g(SH(<p) + a(ip,t),w 0 + t + u(ip,t),t) - 70(7’, 0 - 'yi(<p,t)u(<p,t). 


(2.5) 


The logic in the rewriting (2.3) is that the (linear) approximate dynamics obtained by neglecting 
G is dissipative, with contraction rate proportional to e, thanks to the second condition in (1.4): 
this will allow us to control the full dynamics as a perturbation of the approximate one. The 
approximation obtained by neglecting G is the simplest one displaying dissipation. In principle we 
could have expanded the dynamics at first order both in u and in Sh + a, but the result would be 
qualitatively the same. We now set 


T (tp,t,r):=J dr^i^r') 

and fix wq such that T(<£>, 2n, 0) = T < 0; then we obtain 

u(<p,t) = e er(lM0) [/(</>) + £(tp,t), £(tp,t):=e [ dre er( ^’ T) (jo(<P,t) + G(ip,r)). 

Jo 

The equation for x, if expressed in terms of a, gives, after integration, 

a(<p,t) = £ d Tf(SH(<p) + a((p, t),wq + t + u(lp,t),t). 

Jo 

For t = 2-7T, these give 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


f2 IT 


U(S<p) = e er U(ip) + e / dr e er(¥5,27r ’ r) (7o(y?, r) + G(^,t)) , 


"27T 


H(Sip) = SH(tp) + e / dr f(SH(<p) + a(<p,T),wo +t + u(<p,t),t). 


(2.9a) 

(2.9b) 


It is useful to introduce an auxiliary parameter fi. to be eventually set equal to e, and rewrite 


(2.7) and (2.9a) as 

u(<p,t) = ^ r ^’ t ’ 0) U(<p) + £(<p,t) , £(<p,t) = E [ dr e /xr (' p,t,r ) (yo(<p, r) + G(tp, r)), (2.10a) 


f2 IT 


U(Scp) = Q^ v U{tp) + £ / dr e / " r ^’ 27r,T ^ (70 (ip,r) + G((p,r)), 


(2.10b) 


The idea is to first consider g, as a parameter independent of e, then write the solution in the form 
of a power series in e. with coefficients depending on fi. and hnally show that the (/^-dependent) 
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radius of convergence of the series in e behaves like /x“, a < 1, at small /x: this implies that we will 
be able to take // = e without spoiling the summability of the series. 


Summarizing, we will look for a solution of (1.3) in the form 


with 


x(t) = Stp + Sh(ip) + a(tp,t), w(t) = wo + t + 


£(v,t) = £ [ dr e # * r(v ’ ,t,T) (70(92, r) + G{ip, r)), 

Jo 

r2n 

U(Stp) = e^U (93) + e / dr e ^ r ^’ 27r ’ T ) (70(93, r) + G(ip, r)), 

Jo 

a(ip,t) = s d TF(ip,r), 

Jo 

r2n 

h(Sip) = Sh(<p) + e / d TF(ip,r). 

Jo 


( 2 . 11 ) 

(2.12a) 

(2.12b) 

(2.12c) 

(2.12d) 


and: F(ip,t) := f(x(t),w(t),t), with x(t) and w(t) as in (2.11); G is defined in (2.5); fi must be 
eventually set equal to e. 

2.2 Recursive equations 


The solution to (2.11)-(2.12) is looked for in the form of a power series expansion in e (at fixed fi, 


in the sense explained after (2.10)). Therefore, we write 


1/(93) = fVt/(">( 93), Z(<p,t) = J2e n ^\<p,t), 

n —1 n —1 

00 00 

h(<p) = J2 £n h {n) (<P), a(<P,t) = ^W"^,*), 


(2.13a) 

(2.13b) 


n =1 


n= 1 


and insert these expansions into (2.12). By the analyticity assumption on / and g , we may expand 


(defining q = (q + , q_, q 3) and No = Z D [0, +00)) 

G(<p, t) = Y 0 (A+Mv?) + a+((p, t)) g+ (A -h- (93) + o_(93, t)) 9_ (e^’^ 1/(93) + £(93, t)) 93 , 

qSN 3 

F(<P, 0 = Y F *) (A+MtO + a +( 7 b *)) 9+ (A-MtO + a-(^, i)) Ci ~ (e^ r(¥, ’*’ 0) 1/(93) + £(97i))®\ 
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where * in the first sum denotes the constraint q 7^ (0,0,0), (0,0,1), and h± = h ■ x± and similarly 
for a± (recall that \± and x± are the eigenvalues and eigenvectors of S'); here and henceforth we 
are denoting by • the standard scalar product in M 2 . Moreover 

Gq(<p, t ) := —d q + d q _~d^g(S(p, w 0 + t, t ), F q (ip, t ) := —d q + d q _~ dfif(S<p, w 0 + t, t ), 

q! q! 

where q! := q^lq-lq^l and d± := x± -d x . For future reference, we note since now that the analyticity 
of / and g yields, by the Cauchy inequality, 

IG'qUFql^ 1 , (2.14) 


for some constant Co > 0, uniformly in ip, t (here |q| = q + + q~ + <73). 

Define Z + := No, := Z\No, and F ac{ := x a ■ F q , with a = ±). Setting A := e^ r and plugging 
(2.13) into (2.12), we find for n = 1 

dre /ir(¥, ’ f ’ T) 7o(^,r), 

Jo 

U^\p) = Y A” m_1 / dre^ sm ^’ T )7o(S , >,r), 


m£Z_ 

rt 


a aH < P,0= I dri^o^T), 


p2tt 


h£\<p) = -a Y A a m_1 / drF a ^S m p,r). 


mGZ Q 


Introduce the notation 


/ q + 


Jq,ri(^,T) := JJ(A+/if i:) (^) + aj^for)) 


\i=l 

/ q++q- N 

f II (A-^M + a^^r)) 

\j=<j++i / 


(2.15a) 

(2.15b) 

(2.15c) 

(2.15d) 


(2.16) 


' |q| \ 

H + ^ ni \ip, T)) , 

U=q + +g_+l / 


where n := (m,..., ni q |). Here and henceforth, if q + = 0, the product nSi( A +^+^(¥ , )+ a + r )) 
should be interpreted as 1, and similarly for the other products in the case that q~ = 0 and/or q% = 0. 
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Then, defining Ni :=Zfl (0, +oo), we find, for n > 2, 


II 

Uy> 

[ dre^^’ 
Jo 




(2.17a) 



qeNg 


neN( qi : 

1 77, | —77,— X 





r2n 





U^ n \ip) = 

E *“ m “ 1 j 

1 dr i 

e^r (r Vl : 


) E P^n(S m V,T), 

(2.17b) 


raEZ_ 

0 


qeNg 

neN( q| : 







|n|=n—1 


< 4"W) = 

[ dT £ ^ 

(,q i^Pt T 

) E 

Pq,n((Pi t')i 


(2.17c) 


qe^o 


neN( q| : 

17l| —7T.— 1 




h£Hv) = 

-a E V”' 

r2n 


Fa,q(S m <p,T) E 



- 1 / 

dr E 

P q ,„(S”V,T), 

(2.17d) 


mEZ a 

JO 

qSNg 

neN( q| : 




where |n| := m + ... + ri| q |. 


2.3 Tree expansion and convergence 


We now want to bound the generic term in the series originating from the recursive equations 


(2.17); the goal is to show that the n-th order is bounded proportionally to (Cp, ") , with C > 0 


and 0 < a < 1. We find convenient to represent graphically the coefficients in (2.13) in terms of 


rooted trees (or simply trees, in the following) as in |9j. We refer to [9j Section V] for the definition 
of trees and notations. With respect to the trees in [9] in the present case there are four types of 
nodes. We use the symbols •, O, ■ and □, calling them nodes of type 0, 1, 2 and 3, respectively: 
they correspond to contributions to £, U, a, h, respectively. The constraint * forbids a node of type 
0 or 1 to be immediately preceded by exactly one node of these two types. 

Recall that a tree is a partially ordered set of nodes and lines; the partial ordering relation is 
denoted by y and each line will be drawn as an arrow pointing from the node it exits to the node 
it enters. We call V ( 9 ) the set of nodes and L(9 ) the set of lines of the tree 9. As in [2] we denote 
by vq the node such that vq y v for any node v E V(9): vq will be called the special node and the 
line exiting vq will be called the root line. The root line can be imagined to enter a further point, 
called the root, which, however, is not counted as a node. 

With each node we associate a label r) v E {0,1,2,3} to denote its type and a time variable 
t v E (0, 27t]. With the nodes of types 2 and 3 we also associate a label a v E {±}. Denoting by 
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V(v,w) the path of lines connecting w to v -< w, with both v and w included, set 


m(v) = ^2 mw > 

w&V{v,vq) 


where m v € Z_ if rj v = 1, m v G Z Qti if r/ v = 3, and m v = 0 otherwise. Given a node v, we 
denote by v' the unique node immediately following it; moreover, we let si be the number of nodes 
of type 7j immediately preceding it; if rj = 2,3, we also define sl’ a , a = =L, to be the number of 
nodes w of type r] with a w = a immediately preceding it (i.e. such that w' = v). Finally, we let 
q t , = (</„,+ , q v ,-,q v , 3), with q V; ± = sl }± + s^ ,=t , and q v ^ = + s*. A node v is called an end-node if 

q. ( , = (0, 0, 0), while it is called an internal node if it is not an end-node. 

The node factor A v = A v (ip, t v >, t v ) is defined as 


A, ■= < 


^-m v -l e /ir(5 m W^,27r,r„)^ T v ) e^ s i r ( S ' m(, ’V>'U>, 0 ) ^lv,+ yQv, 

F av , q„ T V ) e A*«ir(S m (’’V,r«,0) A 5«,+ A 9»,- ? 

-a v \-™ v ~ l F av ^{S m{v) ^ r v ) 


H—,v 

1 


if rj v = 0, 

if 7j v = 1, 

if r) v = 2, 
if r/,, = 3, 


where t v i is interpreted as equal to t, while the node integral is 



if rj v = 0, 2, 
if r] v = 1,3. 


With the definitions above, we denote by 0*^ the set of labelled trees with n nodes, rj VQ = 77, and 
the constraint that nodes of type 0 or 1 cannot be immediately preceded by exactly one node of 
type 0 or 1; if 77 = 2, 3, we also denote by 0* ;r?a the subset of by 0* )r) with a VQ = a. Then, one can 
prove by induction that 


22 Val w> u^( (p )= J2 Val ( 0 )’ 

flee;,! 

h£X<P) = £ Val W- 


see; 




(<p,t)= 22 Val ^)’ 


eee; 


eee* 


(2.18a) 

(2.18b) 
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where 


Val(0)= H I V A V , 

vev(8) 

with the integrals to be performed by following the tree ordering, i.e. by starting from the end-nodes 
and by moving towards the root. 


In Figure 2.1 the first order contributions are graphically represented, while the contributions 
of order n = 2 are shown in Figure [2~2| For each node v, the label a v is drawn superimposed on the 
line exiting v, for clarity purposes, while the label m v (to be summed over) is not explicitly shown. 


t (1) (<P,t) = 


VO 


uW(<p)=- 


-o 


VQ 


ai 1} (tp, t) = 


VQ 


ha\<p) = ■ 


-□ 


i;o 


Figure 2.1: Graphical representation of U^(ip), a,a\tp,t) and ha\p)- 


£ (2) M=- 


^ VI 


VO 


U^(p) = 


^ vi 


O 


VQ 


VI 


aa ) (ip) = ■ 

+ 


^0 


-n f 


"VQ ^ Vl 


+ 

+ 


VQ 


VQ 


Vl 


ha' 1 (p) = 


-O 


VQ ^ Vl 


+ 


-Of. 


VQ 


+ 


-O 


H~o ^ v\ 


+ 


-D 




Vl 


Figure 2.2: Graphical representation of t), U^>(p,t), a„\p,t) and ha\ip,t). 


For n = 3 a few contributions to U^(p) are shown in Figure 2.3 All the other contributions 
are obtained by replacing the nodes preceding the special node uo by nodes of a different type, with 
the constraint that if only one line enters vq then it exits a node of type 2 or 3; note that in [9j the 
linear trees at the bottom of the figure were not possible. 
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Figure 2.3: Graphical representation of some contributions to U^ 3 \ip,t). 


Given 9 E 0^, we lot T{9) be the family of labelled trees differing from 9 just by the choice of 
the labels {m v }. Then, using (2.14) and proceeding as in [9} Section VI], we obtain 


£ |Val(0)|<Cr/x-^W, 

6>eT(0) 


(2.19) 


where C\ > 0 is a suitable constant and N*\0) is the number of internal nodes of type in 0; see 
[9] for details. We are then left with bounding Nj (6). 


Lemma 2.1. For all 0 E 0* v one has (9) < 


2 n- 1 


The proof is given in Appendix |a| where it is also shown that such an upper bound on A r ? ! (0) 
is optimal, i.e. there are trees that saturate the inequality. Combining Lemma 2.1 with ( 2.19[ ) and 
recalling that the number of distinct families T{9) in 0* is bounded by C%, for a suitable C*2, we 
get 


max 


{ie (n) (^t)U^ (n) (^)l,k n) (^i)U^ n) Ml} < (2.20) 


where C3 is a positive constant and [•] is the integer part. Eq.(2.20) implies that the radius of 
convergence of the series (2.13) is bounded by £0 = Cj" 1 // 2 / 3 . Therefore, we can take // = e. Holder- 
continuity of U and H can be proved mutatis mutandis like in j§]. This completes the proof of 
Theorem o 


12 




















3 Lyapunov Exponents 

3.1 Conjugation of the tangent dynamics 


In order to compute the Lyapunov exponents, we need to understand how the vectors on the 
tangent space evolve under the interacting dynamics. To this purpose, we set x := {xi,X 2 ,x%) = 


(xi ,x 2 ,w) = (x,w) E T 3 , rewrite ([L3]) as x = f 0 {x,t) + ef{x,t), with / = (/n/2,/3) := (/i,/2,s), 
and write the dynamics on the tangent space as follows: 


y = A t (Sl{x))y, A t {x) = 5{t) [P+ log A+ + P_ log A_] + edf{x,t) , 


(3.1) 


where y E M 3 , S f e (x) = (x(t),w(t)) is the solution (2.11) to (1.3) found in Section [2] P± are the 
projections into the unperturbed eigendirections (x±,0), X± are the corresponding unperturbed 
eigenvalues, and df is the Jacobian matrix of f. 

Integration of (3.1) gives the tangent map DSl : T x T 3 = M 3 —> T s t^ T 3 = M 3 . We denote 


by DS t e (x)y the solution to (3.1) with initial condition y at t = 0 (started at S®(x) = x). For 
t E (0, 27t] we obtain 


ft 


(2 

1 

A 

/ dr df{S T £ {x),T)DS T e {x)y , 

A = 

1 

1 

0 

Jo 


v> 

0 

1 


(3.2) 


We look for a conjugation % : T 2 x M 3 —> T 3 x M 3 , 

U : {ip, Y) —► (H{ip), (1 + K{ip))Y), H{ip) = W{ip)), 

where 1 is the identity in M 3 and K(ip) is a 3 x 3 matrix, such that, by setting 

*!(*, y) = (SUx), DSi(x)y), S 0jE (<p , V) = {Sip, (A + N{ip)) Y ) 
for a suitable 3x3 matrix N{<p), to be determined, one has o % = T~L o Sq e , that is 
^ n (H{ip), (1 + K{ip)) Y) = {H{Stp), (1 + K{Sip)) (A + N{<p)) Y). 


(3.3) 


(3.4) 


Of course only the part of (3.4) involving the tangent dynamics has still to be solved, so we study 
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the conjugation equation 


DS 2 /(H(<p)) (1 + K{<p)) Y = (l + K(Sip)) (A + N(<p)) Y. 


(3.5) 


The matrix N(ip) will be taken to be diagonal in the basis |y + , y_, r/ 3 }, where y± = (x±,0), 
and y 3 = (0,0,1). Then in the basis {y + ,y_, y 3 } one has 


K(<p) = 


while the matrix A takes diagonal form, with values A = (A+, A_, 1) := (A+, A_, A 3 ) along the main 
diagonal. From now on we shall use this basis, and implicitly assume that the indices run 

over the values , 3, unless stated otherwise. 


A++M 

k+-(<p) 

fc+3(^)\ 



0 

0 ) 

k~+{<p) 

k—(<p) 

k-3(<P) > 

N(<p) = 

0 

v -(<p) 

0 

\fa+(<P) 

fa -(¥>) 

fa3 (<p)J 


l o 

0 

^3 (</>)/ 


By setting y^) : = (1 + K((p))y i , we obtain from (3.5) 


DS? r (H(<p),W(<p))y i ((p) = A i(<p)yi(S(p), A i(<p) := A i + i/ifa). 


(3.6) 


Note that, given a solution yi(<p) of (3.6), then also y'i(<p) := h(<p)yi(tp) is a solution with A i(ip) 


replaced with A'((/?) = ( k(tp)/li(Sip )) A i(tp), where k((p) are non-zero functions from T 2 to M. There¬ 
fore, with no loss of generality, we can require the diagonal elements of K{<p) to vanish: hence K(<p) 
will be looked for as an off-diagonal matrix. 


3.2 Recursive equations 


We look for the conjugation in the form (3.3)-(3.4), with 


kj(<p) = M<P) = 5^e n ^ n) (^). 


(")/ 


71=1 


72—1 


Equations (3.2) and (3.5) give 


r2n 


A (1 + K(ip)) + e / &Tdf(Sl{H{v)),T)DSl{H{v)){t + K{v)) 


(3.7) 


= A + N{(p) + K(S<p) A + K(S<p)N((p), 


(3.8) 
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where DSg(H(<p)) can be computed iteratively via (3.2) as 


n> 1 J ° J ° 

• f dT n df(S£ n (H(tp)),T n )A:=(t + Wl(tp,T))A. 


(3.9) 


At first order in e, (3.8) gives (with i, j £ {±,3}) 

r2n 

~ ^j k ij\Sip) + Aj J dTdjfi(Stp,w 0 + T,T) = v^\ip)Sij. (3.10) 

In particular, setting i = j and recalling that fcW is off-diagonal, we find 


f'2n r2n 

t '±\<p) = ^± dTd±f±(Sip,w 0 + t,t), vf\(p)= dTd w g(S(f,w 0 + T,r) = T, 

1 0 ./o 

.( 1 )/ 


while, for z 7^ j, we have to solve recursively (3.10) for fcb (</?), the result being (if a = ±): 


&£-«(¥>) = ~« 2 


m£Z 0 


A-r 


—m—1 


dTd- a f a (S m+1 ip,w 0 + T,T), (3.11a) 


f27T 


k aU‘f) = ~ a A “ m 1 / d^/aOS^Vwo + r,r), 
mez a 

r2ir 

4!-a(¥>) = -a X] A -« 1 / dr9_ Q / 3 (5' m+ V,w;o + r,r). 


(3.11b) 

(3.11c) 


mGZ 0 


In order to compute the higher orders, we insert (3.9) in the left side of (3.8), thus getting 


AK(ip) - K(Stp) A - N(<p) - K(S<p)N(ip) + 97t(<^)A(l + K(<p)) = 0, (3.12) 


where 9JT(<^) := Wl((p,27r). Note that, according to (3.9), VJt(ip) is expressed as a series of iterated 
integrals of df (Sg(H(ip )), r), where df is analytic in its argument. Therefore, the power series 
expansion in e of 9Jl(</?) can be obtained (and its n-th order coefficient can be bounded) by us¬ 
ing the corresponding expansions for the components of Sg(H(ip)) = ( SH(ip ) + a((p,r),wo + r + 
e Mr(¥>,T,o + €{tp,t))\ here the functions H(<p), a(ip,r), U(cp) and £(<£>, r) are as in (| 2.13 [) with 


15 










the coefficients given by (2.18) and bounded as in (2.20). We write 


91 l(<p) = '^E n Vn (n) (p). 


(3.13) 


n— 1 


By using the very definition of 9JT(y?) and the bounds (2.20), it is straightforward to prove that 

||91t (n) (^)|| := max |^j(¥>)| < Q/x -[(2n “ 1)/31 (3.14) 


for a suitable C4 > 0. Now, if n > 2, the diagonal part of (3.12) gives 




>M 


On)/ \\ 1,(^2)/ 


ni+n2=n j —±,3 

ni,7i2>l 


while the off-diagonal part can be solved in a way similar to (3.11), i.e., if a = =L, 

—m 


= -&K 1 Y 


mG Z Q 


A, 


A-c 


M { Za(S m <p)\-a + QZa(S m <p) 


(n) 


ni-\-ri2=n 

ni,n 2 >l 


l'=±,3 


In the simple case that f = (0,0, g), v+(<-p) = v-(<p ) = 0, while 


/•2n rri rr n - 1 2 -e 

1 + I*^) = 1 + V e n / dn dT2 ■■ dT f[ d w g(S<p, w(Ti),n) = e £ f° «^(S^(r),N 

n> 1 J 0 J 0 J 0 i =1 

thus recovering the formula for ^4(27 t) = log(l + 1/3 (<£>)) given in |9} Section VII], 


(3.15) 


(3.16a) 

(3.16b) 

(3.16c) 


k^) = -a\~ 1 Y A+ 

mGZ Q 

=-«E *-» K‘b(s"V)A-« + ei;l(s m V ) 

mGZ Q 

where we have set 

sg>(9):= E E (M7) 
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3.3 Tree expansion and convergence 


In Figure 3.1 and 3.2 we give a graphical representation of (3.15) and (3.16a), respectively. The 


representation of (3.16b) and (3.16c) is the same as in Figure 3.2, simply with the labels (a,—a) 
replaced by (a, 3) and (3, a), respectively. 



Figure 3.1: Graphical representation of u- n \(p). 




(m) 


(”2) 

3.2: 


Graphical representation of 


C-V)- 


To iterate the graphical construction and provide a tree representation for both N(tp) and 
we need a few more definitions. We identify three types of principal nodes, that we call of type K, 
N and 9J1, and represent graphically, respectively, by □ , ■ and 9 . With any such node v, we 
associate a label rj v E {K,N,Wl} , to denote its type, and two labels i v ,j v E {+,—,3}, which will 
be drawn superimposed to the line exiting v\ if v is of type N, then i v = j v , while if v is of type K, 
then i v ^ j v . A node is of type 911 if and only if it is an end-node. Furthermore, with each node v 
with rj v = 9J1, we associate a label n v E Ni, while we set n v = 0 for all nodes v with r] v ^ SOI; with 
each node v with rj v = K, we associate a label a v E {+, — } such that either i v = a v or j v = —a v 
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(recall that if v is a nodes of type K then i v ^ j v , so that either i v or j v are = ±), and a label 
m v E Z a „; if v is of type N or 911 we define m v = 0. 

If denotes the number of lines entering v and si the number of lines of type r] entering v, we 
have the constraints s v E {0,1, 2} and si E {0,1}. Moreover: s v = 1 =>■ s = 1; s v = 2 =>■ = 1. 

If = 1 and is the node immediately preceding v on has i v = i Vl and j v = j Vl . If s v = 2, let 
vi,i >2 be the two nodes immediately preceding v, if s = 1, with no loss of generality we assume 
that v\ is of type 911 (so that V 2 is of type K); if = 1, with no loss of generality we assume that 
v\ is of type K (so that V 2 is of type IV); in both cases we impose the constraints that i v = i Vl , 
jv = jv 2 and j Vl = i V2 . Denoting by v' the node immediately following v, we set q v '■= s^,, and 

P(v) = ^2 ( m w + qw)- 

w&V( v,vo) 

We are finally ready to define the node factors associated with the nodes: 


Ay (<^) 


1, if rj v = N, 

< -a v Xr\x iv /\ jv )-"*(-i)‘Z, if Vv = K, 

x jv an<zl(SPM<p), if r h = fm. 


(3.18) 


Then, by iterating the graphical representation in Figures 3.1 and 3.2 we end up with trees like 


that in Figure |3l3| for i j: note that the end-nodes are all of type 91T. If i = j the only difference 
is that the special node vq is of type N. 



Figure 3.3: A tree contributing to with n = 11. 
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With the definitions above, we denote by 0** ); ^ ^ the set of labelled trees such that Y!veN(6) n v = 

n, r] VQ = rj. and the constraints and properties described above. Then it is straightforward to prove 
by induction that 




(<p)= E Vai w, k ij( ( p)= E Val ^)> 


(3.19) 


6»e0; 


n,N,(i,i ) 


oee 


n,K,(i,j ) 


where 


Val(0):= H A v . 

vev(9) 


(3.20) 


Given 6 G 0** ?; ^ we let T{6) be the family of labelled trees differing from 8 just by the choice of 
the labels {m„}. Then, using (3.14), it is easy to see that 


E |Val(0)|<( T 

0eT(0) 


A + \JVfW 


- A 


-2 


c m U2n-E(0))/3] 


which immediately implies that 


kg’M . 4"V) < c 5 ”m l(2 ”- 1)/3 ', 


for a suitable constant C 3 > 0. Therefore, the radius of convergence in £ of the series for K(tp) and 
N e ((p) is proportional to ^r~ 2 / 3 , which allows us to fix eventually e = fi. 

The Lyapunov exponents A* are the time average of the quantities log Xi(ip). However, if S’f 7I jfi 
denotes the restriction of S 2n on the attractor H, the dynamical system (H, £'| 7r |H) is conjugated to 
an Asonov system and hence it is ergodic: therefore time-averaged observables are (^-independent. 
Furthermore, there exists a unique SRB measure fj, £ such that 

A± = J d ii e (tp) log (A± + v± (<p)) , A° = j d n e {ip) log (1 + v 3 ((p)) . (3.21) 


The measure fi £ can be computed by reasoning as in j8j Chapter 10]. Let be a Markov partition 
for S on T 2 and set W(^Po))- Call X e the symbolic code induced by the Markov 

partition and denote by a_ the symbolic representation of a point x£(l, i.e. x = X £ (a). Then 
the expansion rate of S' 2?r |H along the unstable manifold of X £ (a) is e A+ (-\ where 


A + (a) = log (A+ + v + (X 0 (a))) 


|w + (X 0 (rg))| 
|w + (X 0 (u))| ’ 
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with t denoting the shift map and w_|_(<£>) = (1 + K((p))y+. If p, £ denotes the Gibbs distribution 
for the energy function A + (q_) (see [8J Chapter 5]), then the SRB distribution hsrb for the system 
(f2, Sf 7r |fl) is fi e and can be computed accordingly (see 0 Chapter 6]). 

Moreover, by construction, A+(a) is analytic in e and Holder-continuous in a. Therefore, for 
any Holder-continuous function F : T 2 —> M, the expectation value 

Hs{F) := J d /j, e (ip)F(ip) 


is Holder-continuous in e for 0 < e < £o- If F is analytic in e, then there exists a positive constant 
ep, depending on F. such that /i e (T) is analytic for |e| < ep < eo- In particular the Lyapunov 
exponents (3.21) are analytic in e and, from (3.21), one finds 


A 0 = / d He(tp) log (1 + v 3 (tp)) = j d fjL e (y) (eT + O (e 2 )) = eT + O (e 2 ) . 


This completes the proof of Theorem 


1.2 


4 Perspectives: transition from fractal to smooth(er) synchroniz¬ 
ation 

In this section, we discuss informally some of the consequences of our main theorem, and formulate 
a conjecture about the transition from fractal to smooth(er) behavior, which is suggested by our 
result. From Theorem o we know that the surface Q of the attractor is Holder continuous, but 
we do not have any control on its possible differentiability. This means that our attractor may be 
fractal, and we actually expect this to be the case for e positive and small enough. An analytic 
estimate of the fractal dimension of the attractor in terms of the Lyapunov exponents is provided 
by the Lyapunov dimension Dl, which is defined as follows. Consider an ergodic dynamical system 
admitting an SRB measure on its attractor, and let Ai > A2 > ■ ■ ■ > A n be its Lyapunov exponents, 
counted with their multiplicities. Then, 

k 

= t + (4 ' 1) 
k 

where 1 < k < n is the largest integer such that ^A i > 0. The Kaplan-Yorke conjecture [6| 

2 — 1 

nn states that Dj j coincides with the Hausdorff dimension of the attractor (also known as the 
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information dimension , see, e.g., [3 Chapt.5.5.3] for a precise definition). In this section, we take 
Dl as a heuristic estimate of the fractal dimension of the attractor, without worrying about the 
possible validity of the conjecture (which has been rigorously proven only some special cases, see 
e.g. |5]). 

Specializing the expression of Dl to our context, we find that, for e sufficiently small, 

Dl = 2 H j- : = 3 + e- - b R 2 (e), (4.2) 

I A-1 log A+ 


where i?2(^) = 0(e 2 ) is the Taylor remainder of order 2 in e, which is computable explicitly in terms 
of the convergent expansion derived in the previous sections. Note that T < 0, so that Dl = Dl{e) 
is smaller than 3 (as desired) and is decreasing in e, for e small. Therefore, combined with the 
Kalpan-Yorke conjecture, ( |4.2| ) suggests that the attractor is fractal for e small, and its fractal 
dimension decreases (as expected) by increasing the strength e of the dissipative interaction. 

It is now tempting to extrapolate (4.2) to larger values of e (possibly beyond the range of validity 
of Theorem 1.2), up to the point where, possibly, the relative ordering of Ao and A + changes. In 
the simple case that / = 0 (which is the case considered in (9j), the Lyapunov exponents A± are 
independent of e: A+ = —A_ = logA+. Therefore, on the basis of (4.2), we conjecture that by 
increasing e the Hausdorff dimension of the attractor decreases from 3 to 2 until e reaches the critical 
value e c , where Ao(e)| £=£ = A_. Formally, this critical point is e c = (log A+)/(—T)+(higher orders), 
the higher orders being computable via the expansion described in the previous sections. For e > e c , 
we expect the attractor to be a smooth manifold of dimension two. The transition is illustrated in 
Fig |4.1| and 4.2 for the simple case that / = 0 and g(ip, w, t ) = sin(u; — t ) +sin((^2 + w + t), in which 
case the expected critical point is e c ~ 0.153. 




Figure 4.1: / = 0 and g(ip, w, t ) = sin(w; — t) + sin(<^2 + w + t). a): e = 0.05. b): e = 0.1. 

If / 7^ 0, on the basis of numerical simulations, the attractor does not seem to display a transition 
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from a fractal set to a smooth manifold. Still, for suitable choices of / ^ 0, we expect the attractor 
to display a “first order phase transition”, located at the value of e where Ao(e) = A_(e), to be 
called again e c . At £ = e c , the derivative of the Hausdorff dimension of the attractor with respect to 
e is expected to have a jump. A possible scenario is that the attractor is fractal both for e < e c and 
for e > £ c , but it is “smoother” at larger values of e, in the sense that its closure may be a regular, 
smooth, manifold of dimension two. An illustration of this “smoothing" mechanism is in Fig|4.3[ 



It would be interesting to investigate the nature of this transition in a more quantitative way, 
by comparing a numerical construction of the attractor with the theory proposed here, obtained by 
extrapolating the convergent expansion described in this paper to intermediate values of e. Such a 
comparison goes beyond the purpose of this paper, and we postpone the discussion of this issue to 
future research. 
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Figure A.l: A tree of order n = 11 for which the bound on Ac 1 (0) is optimal. 


A Estimate of Nl(0) 


One has Nl(9) = 0 for n = 1 and A r l (0) < 1 for n = 2. Given a tree of order n > 3 one proceed by 
induction. Assume that All (O') < (2 n' — l)/3 for all the trees 9' of order n' < n — 1 and consider 
a tree 9 of order n. Let no be the special node of 9, and call 6i,..., 6 S the subtrees entering no, 
with s > 1. If no is a node of type 1, then N^(6) = 1 + Ylk=i A;* (9k) < 1 + Efc=i(^ n fc — l)/3 < 
1 + (2 (n — 1) — s)/3 = (2 n + 1 — s)/3, so that the bound follows for s > 2. If s = 1, then the node 
preceding no cannot be of type 1. Call 9[,... ,9' , the subtrees entering m, with s' > 1. Then one 
has Nl(9) = 1 + J2(=i N l(6'k) < 1 + ( 2 ( n ~ 2 ) - s')/3 = (2n - 1 - s')/3 < (2n - l)/3 for all s' > 1. 
Finally if no is not a node of type 1, then it has not to be counted and the argument follows by 
using the inductive bounds for the subtrees entering no- 

Moreover the bound in Lemma |2.1| is optimal. Indeed there are trees 9 of order n such that 
A'* 1 (9) = (2 n — l)/3. Define recursively the level £(v) of a node n by setting £{v) = 0 if n is an 
end-node and i{y) = p if at least one line entering v exits a node w with level £(w) = p — 1. Then 
consider a tree in which all nodes except the end-nodes are circles (thais is of type 1 or 2) and have 
two entering lines except those with level 1 which have only one entering line; see figure |A.l for 
an example with n = 11 (J5] means we can have any kind of square node). For such trees one has 
n = 2 k — 1 + 2 k ~ 1 , where k = £(vq), 2 k — 1 is the number of internal nodes and 2 fc_1 is the number 
of end-nodes. Hence (9) = 2 k — 1 = (2 n — l)/3. 
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